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Abstract: It has been realized for a long time that knowing the ij and r]' wave functions 
in terms of quark and gluon components probes our understanding of nonperturbative 
QCD dynamics. Great effort has been given to this challenge, yet no clear picture has 
emerged even with the most recent KLOE data. We point out which measurements would 
be most helpful in arriving at a more definite conclusion. A better knowledge of these wave 
functions will significantly help to disentangle the weight of different decay subprocesses in 
semileptonic decays of D~^, and mesons. The resulting insights will be instrumental 
in treating even nonleptonic B transitions involving r] and rj' and their CP asymmetries; 
thus they can sharpen the case for or against new physics intervening there. 
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1 Introduction 

The question of r] — ij' mixing ^, i.e., how their wave functions are composed of SU{3)fi 
singlet and octet qq components, goes back to the beginning of the quark model era [1-9]. 
With the advent of QCD it became even more involved, since QCD brought with it more 
dynamical degrees of freedom, namely, gluons, which can form a second class of SU{3) fi 
singlets. Determining t] — rj' mixing is thus an intriguing element in understanding QCD's 
nonperturbative dynamics. Lattice QCD's attempts to establish theoretical control over 
this mixing are still in their infancy [10, 11]. Showing that there is a purely gluonic compo- 
nent in the r] and/or rj' wave functions would establish for the first time that gluons, which 
have been introduced to mediate the strong interactions and whose presence as indepen- 
dent degrees of freedom has been demonstrated as progenitors of jets in 'hard' collisions, 
play an independent role also in hadronic spectroscopy. In section 2 we introduce basic 
notions relevant for rj — rj' mixing, while in section 3 we review the somewhat ambivalent 
findings from several phenomenological studies. Armed with this knowledge we discuss 
weak D and B decays producing rj and rj' mesons in section 4 and what the observed rates 
can tell us about the underlying quark level transitions; we comment briefly on how the 
structure of the rj and rj' wave functions affect CP asymmetries in the channels B^ — )• rj'Ks 
and Bd — >• r]Ks- Finally in section 5 we present a summary and outlook. 

^The term "mixing" is often used when oscillations, e.g., B° — B'^ are involved; however with oscillations 
one has a nontrivial time evolution, but not for rj and ij' mixing. 
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2 T] — 7]' Mixing 



Based on approximate QCD flavor SU{2)) fi symmetry, the mixing of the rj and r]' mesons 
can be described in two different bases: 

1. the SU{3)fi singlet and octet components jryo) = -^\uu + dd+ss) and \7js) = -^\uu + 
dd — 2ss), respectively: 




cos Op — sinOp 
s'mOp cos 9p 




(2.1) 



the quark-flavor basis with |r/g) = + dd) and |r/s) = \ss) 




cos <pp —smcpp 
sin0p cos (j)p 




(2.2) 



As long as state mixing is regarded, one may freely transform from one basis to the other; 
the two parametrizations are related through 

6»p = 0p - arctan\/2~ (/)p - 54.7° (2.3) 

In the SU{3)fi symmetry limit, 6p = 0, and (pp takes the so-called 'ideal' value (j)p = 
arctan\/2 ~ 54.7°. 

Just for orientation: the quadratic [linear] Gell-Mann Okubo (GMO) mass formula 
points to Op ~ -10°, 0p ~ 44.7° [Op ~ -23°, <j)p ~ 31.7°]. 

The mixing schemes have been analyzed in the context of chiral perturbation theory. 
On lattice, it is not an easy task to study r] and r]' , as experienced in the last decade of 
attempts. The RBC-UKQCD Collaboration has reported a pioneering calculation of the r] 
and ij' masses and mixing angle of Op = —14.1(2.8)° using Nf = 2 + 1 flavor domain wall 
ensembles on an Iwasaki gauge action [10]. Their results show small octet-singlet mixing, 
consistent with the quadratic GMO within the large statistical errors. Masses and mixing 
angle of the rj and rj' have also been calculated by the Hadron Spectrum Collaboration [11], 
using lattice QCD with unphysically heavy light (up, down) quarks and a single lattice 
spacing: their estimate value is (j)p = 42(1)°. The large value of the mixing angle (pp in 
the pseudoscalar sector, with respect to other ones (e.g. the vector mesons \uj >~ jry^ > 
and 4> ~ \r]s >, with a mixing angle (pv = (3-4 ± 0.2)° [12]) is expected, because of the 
additional mixing induced by the axial U{1) anomaly ([13] and references therein). 

In the 1990s the possibility of a single angle description being inadequate started to be 
considered. Several papers [13-21], based on theoretical studies as well as on comparison 
with data, pointed out that the pattern of SU{3) fi breaking requires a description in terms 
of two angles. Phenomenological analyses have often involved weak decay constants /^(,), 
defined by the relation < 0|A^|7/^')(p) >= if"'^,)P^- In the octet-singlet basis a = 8,0 

8 

and are the octet and singlet axial-vector currents; in the quark-flavor basis, a = 
q, s and Aju^ are the nonstrange and strange axial-vector currents. Because of SU{3) fi 
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breaking, the mixing of the decay constants does not necessarily follow the same pattern 
as the state mixing (see e.g. [20, 21]). For completeness, we report here the most general 
parametrizations involving two independent axial- vector currents and two different physical 
states: 



/' fv \ = ( /8 cos 08 -/osin^o' 
4' /°' / 1/8 sin 08 /o cos 00 



(2.4) 



2. 




(2.5) 



We observe that in Eq. (2.4) as in Eq. (2.1) the angles are chosen in such a way that 
9p = 9s = Oq = corresponds to the SU{3)fi symmetric world. As before any expression 
in one scheme can be translated into the other one in a straightforward mathematical 
way. However different dynamical implementations of SU{3)fi breaking suggest a different 
ansatz; for example it has been suggested that attributing SU{3)fi breaking to Okubo- 
Zweig-Iizuka (OZI) violating contributions leads to (pg ~ (pg, recovering a description in 
terms effectively of a single angle in the quark-flavor basis [13, 19]. As it is well known, the 
OZI rule leads to a suppression of strong interaction processes where the final states can 
only be reached through quark anti-quark annihilation. In the octet-singlet basis, instead, 
the differences in 9 may be sizable, and most analyses find the range 9s — 9o^ [—19°, —12°] 
([13, 16, 18, 19, 21] and references therein). In this respect, the quark-flavor basis plays a 
privileged role; we will use such a basis in the following, assuming a single mixing angle 
4'P = 4'q = 0s; that correspond to Eq. (2.2). We can see from Eq. (2.5) that under this 
assumption the decay constants follow the same pattern of particle state mixing. 

The plot thickens further still in QCD, for one can form an SU{3)fi singlet not only 
from quark- antiquark combinations, but also from pure gluon configurations with the sim- 
plest one being a gg combination. Since in general all components compatible with the 
quantum numbers of a state can appear in that state's wave function, there is no a priori 
reason why the rj and rj' wave functions could not contain such configurations. On gen- 
eral grounds they will contain also cc (or bb) components, but probably on a significantly 
smaller level, since the mass scale for gluonic excitations is presumably lower than the J/ip 
mass; therefore we will ignore cc (and bb) admixtures in our subsequent analysis. Using 
the quark-flavor basis, we write down [5]: 

^ Xr,\vg)+Yr,\rjs)+Zr,\gg) (2.6) 

One would expect the heavier rj' to contain a higher dose of gluonic components than the 
rj, which is also mainly an SU{3)fi octet. Setting to zero is presumably a pragmatically 
sound approximation. In [22] the authors use a number of parameterization schemes to 
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analyze J/ijj and if)' decays into vector and pseudoscalar mesons; in most cases they find 
a value for the gluonic content of r/ compatible with zero, with an exact numeric value of 
Z'^/Z'^i that is strongly model dependent and ranges from 10-^1 to 0.08. Reference [22] 
also presents a framework, based on old perturbation theory, that allows a much higher 
gluonic content in r/, that is Z'^/Z'^, ~ 1. This result is inconsistent with the analysis of 
the same data made in [23], where = is assumed. 

In the following, we use the approximation = 0, Z^/ ^ and we parameterize the 
two orthonormal states in terms of 0p plus an additional mixing angle (jjc'- 

\vi') ~ cos<^Gsin0p|?7g) + cos(^GCOS0p|r/s) + sva.(j)G\gg) 
It?) ~ cos(i)p\7]q) - sm(j)p\i]s) (2.7) 

As already mentioned it is unlikely that lattice simulations of QCD will determine the 
■q and r}' wave functions in the near future. Phenomenological studies are thus our only 
recourse. Several such analyses have been undertaken recently: while their findings are not 
inconsistent, their messages are ambivalent, as we will discuss in the next section. 

3 Phenomenological Studies oi t] — t]' Mixing 

There are three classes of electromagnetic and strong transitions that can provide infor- 
mation on the mixing angles and the gluonic content: 

• Radiative vector and pseudoscalar meson decays: 

ip' , tp, (j) ^ vs. 'jr] 
p, u ^-yr] 

?/ juj, -fp (3.1) 

• Decays into two photons or production in 77 collisions: 

v' ~^ 77 'v^- f] ^ 11 (3.2) 
11 ^ f] vs. 77 — )■ rj' (3.3) 

• Decays of ip into PV final states with the vector meson acting as a 'fiavor filter': 

ip ^ p/uj/(j) + 7] VS. rf (3.4) 

3.1 Present Status 

Recent papers on the glue content of the rf by KLOE [24] and Li et al. [22] have motivated 
other studies of a range of different processes [23, 25, 26]. Escribano, Nadal [25, 27] and 
Thomas [23] have analyzed all processes of Eq. (3.1), (3.2), (3.3) and (3.4). The old and 
new analysis from KLOE [24, 26] and Escribano and Nadal [25] refer to those processes of 
Eq. (3.1) whose dynamical scale is below 1.02 GeV (that is, including but excluding ^ 
and ■0' decays), while Li et al. [22] have analyzed the ones above 1.02 GeV. In Table 1 we 
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Analysis 


(/)p (Ansatz Z^, = 0) 


KLOE 

Escribano I 
Escribano II 
Thomas I 
Thomas II 

Thomas I with form factors 
Thomas II with form factors 


(41.3±0.3,tat±0.7,j,,)° 
(41.5 ± 1.2)° 
(42.7 ±0.7)° 
(41.3 ±0.8)° 
(41.7 ±0.5)° 
(41.9 ± 1.1)° 
(42.8 ±0.8)° 



Table 1. Fit values for the rj — -q' mixing angle as inferred by different authors from radiative 
decays of vector /pseudoscalar mesons below 1.02 GeV, assuming = 0. Only the KLOE analysis 
includes also constraints from rj' — >■ 77. / labels the results from the analysis without including the 
latest data on (j) ^ rj'^ (KLOE) and (p,w,0) — >■ 777 (SND), while II indicates the same analyses 
performed including them. 



Analysis 


4>p zl 


KLOE 
Escribano I 
Escribano II 
Thomas I 
Thomas II 

Thomas I with form factors 
Thomas II with form factors 


(39.7 ±0.7)° 0.14 ±0.04 
(41.4 ±1.3)° 0.04 ±0.09 
(42.6 ±1.1)° 0.01 ±0.07 
(41.3 ±0.9)° 0.04 ±0.06 
(41.7 ±0.5)° 0.04 ±0.04 
(41.9 ±1.1)° 0.10 ±0.06 
(41.9 ±0.7)° 0.10 ±0.04 



Table 2. Fits allowing for a gluonium component using radiative decays of vector /pseudoscalar 
mesons below L02 GeV. Only KLOE analysis includes also constraints from rj' ^ . I again labels 
the results from analyses without including the latest data on (f> rj'j (KLOE) and (p,LU,(l)) — > 777 
(SND), while II indicates the same analysis performed including them. 

have summarized the results of [23-25], based on radiative decays of vector /pseudoscalar 
mesons below 1.02 GeV. 

The KLOE analysis also includes constraints from tt^/t]' — t- 77, according to the pre- 
scription of Ref. [28]. These results are obtained by including vector-pseudoscalar wave 
function overlaps, assuming the r/(') to be a pure qq state, i.e. = 0, and the dependence 
of the decay widths on the mixing angle as in [25]. 

We see that the different analyses yield very consistent values for the mixing angle, 
namely (pp — 42°, which happens to be close to the value suggested by the quadratic GMO 
mass formula. Including the latest data from KLOE [24] and SND [29] does not cause a 
significant shift. 

In Table 2 results from the same studies are listed, now allowing for a gluonic compo- 
nent in T]' , i.e. Z^, / 0. 

The different analyses again yield consistent values for the mixing angle with (jyp — 42° 
with only KLOE finding a somewhat smaller number. As before the latest data from KLOE 
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and SND do not cause a significant shift. Yet wliile the numbers given for the size of a 
gluonic component are not truly inconsistent considering the stated uncertainties, they 
seem to carry an ambivalent message: while the first and last studies ~ listed as "KLOE" 
and "Thomas with form factors" - point to a significant gluonic component, the others 
do not. We can understand some of the differences. As explained around Eqs. (2.4) and 
(2.5) we think that assuming the mixing of the decay constants to follow the same pattern 
as state mixing is an oversimplification. Only Thomas has gone beyond this assumption, 
and when he includes the form factors he finds some intriguing evidence for a gluonic 
contribution. 

The form factors included in "Thomas" are phenomenological Gaussians, whose aim is 
to introduce a momentum dependence for exclusive processes. In order to understand why 
the findings from "KLOE" and "Escribano/Thomas I-II" for the gluonic content in Table 
2 are as different as they appear (for neither analysis allows for different form factors), 
we can offer one comment, though: only "KLOE" includes r( — t- 77, and that observable 
pushes up the value of Z^/, as pointed out by Thomas. 

In fact, the above theoretical discussion has prompted the KLOE Collaboration to 
perform another fit [26], updated by using the branching ratio values from PDG 2008 
[30], the more recent KLOE results on the w meson [31] and using a larger number of free 
parameters, as suggested by [23, 25]. The fit has been performed in the two cases: imposing 
the gluonium content to be zero, that resulted in (\)p = (41.4 ± 0.5)°, or allowing it free, 
giving (pp = (40.4 it 0.6)°. KLOE new results confirm the gluonium content of rj' at 3a 
level with Z^, = 0.115 it 0.036, in contrast with "Escribano/Thomas I-II" values in Table 
1. Therefore, the actual difference between "Escribano/Thomas I-II" and KLOE values 
appears due to the inclusion in the latter of rj' — )• 77. 

The comparison presented above pointed out that decays into two photons can play a 
key role in the mixing parameters determination. They can be exploited also looking at 
the inverse processes, namely, the production in 77 collisions. 

The L3 Collab. at LEP has published [32] the measurement of the radiative width 
T(r]' — )• 77) produced via the collision of virtual photons, in the reaction e~^e~ — )• e"^e~7*7*, 
7*7* — )• r]' , T]' — )■ TT^TT"^, using data collected at centre-of-mass energies ^/s ~ 91 GeV. 
They compare the photon-meson transition form factor with a model by Anisovich et al. 
[33], that allows a variable admixture of gluonic content, from 0% to 15%. The central 
values of L3 data points favour a low gluonium content, but the whole interval is allowed 
within the large errors. 

Before L3, the same e^e~ — )• e~^e~r]' reaction had been performed at lower energy 
e^e" colliders, by using various rj' decay channels (see Refs. in [32]). Let us review 
some old measurements of the radiative widths T{i]^'^ — t- 77) used to evaluate the mixing 
angles. These estimates did not consider the possibility of gluonic content and refer to the 
octet-singlet basis and the single angle approximation, whose limits have been discussed 
in Sect. 2. To facilitate the comparison, we have quoted the results in the flavor basis, 
using the relation (2.3). The observation of r] meson production from 77 fusion has been 
reported in a 1983 Rapid Communication by the Crystal Ball Collab.; the given mixing 
angle reads (pp = 37.1° it 3.6° [34]. In 1988 they published the radiative widths for vr", r] and 
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rl and determined mixing angles from the experimental averages, finding = 32.3° it 1.2° 
[35]. Two years later both the MD-1 [36] and the ASP Collaborations [37] presented the 
measurement of the rj, r]' — t- 77 widths, with results in agreement within the errors. The 
ASP Collaboration calculated the pseudoscalar mixing angle (j)p = 34.9° ±2.2° [37]. While 
these values are compatible among them, they appear to fall significantly below those in 
Tables 1 and 2. 

A new surge of experimental data and updated analyses is strongly needed. The 
BABAR Collaboration has led the way presenting recent studies on the 77* — )■ rj^'^ transi- 
tion form factors in the momentum transfer range from 4 to 40 GeV^ [38]. They compare 
measured values of the rj^'^ form factors with theoretical predictions and data for the vr'^ 
form factor by using the description of rj — rj' mixing in the quark-flavor basis (2.2). They 
assume no gluonic admixture and a mixing angle (j)p = 41° . The dependence on the trans- 
fer momentum of the form factor for the \ris) state is different from the QCD prediction 
[39] of the asymptotic distribution amplitudes; data points are systematically below the 
theoretical curve. Because of the strong sensitivity of the result for the \r]s) state to mixing 
parameters, an admixture of the two-gluon component in the rj^'^ meson cannot be excluded 
as a possible origin of this discrepancy. 

A new investigation is being performed by KLOE, from the analysis of off-peak data, with 
integrated luminosity L = 240 pb~^, already on tape, devoted to the measurement of the 
77—7-77 rate. The off peak analysis, at ^/s = 1 GeV instead of y/s = 1.02 GeV, allows 
to reduce the main background, coming from resonant contributions (p — t- 777. After the 
full selection, the data set consists of 600 77—7-7/ with 77 —7- vr+vr^vr*^ and 900 77—7-7/ 
with T] —7- vr'^vr'^Tr^; the cross section 17(77 —s- 77) at 1 GeV is under evaluation [40]. The 
upgraded KLOE detector (KLOE-2) will be suited for taking data also at energies away 
from the (j) mass. Taggers designed to detect the outcoming e~^e~ are being inserted into 
the KLOE detector, to provide a better background rejection without going off peak and 
allow precision measurements of the 77 cross section. There is a proposal to increase the 
DA$NE energy up to y/s ~ 2.5 GeV; however, a run at y/s ~ 1.4 GeV is already enough 
to measure the 7/' decay width [41]. 

Starting in September 2009, the Crystal Ball at MAMI has undertaken a huge upgrade, 
with an increase of the MAMI beam energy and the construction and assemblage of a new 
tagging device; one reason of the upgrade is a measure of rj' — >• 77 branching ratio [42]. 

The quoted measurements of the width are obtained with the QED process e"^e^ — )• 
e^e~7*7* — >• e~^e~r]. The 2010 PDG average is taken from such experiments and gives 
T{r] —7- 77) = 0.510 lb 0.026 KeV. The error on the average is 5%, while the errors in 
individual experiments range from 8% to 25%. There is a different type of measurement 
of T{r] —7- 77), not included in the 2010 PDG average, based on the Primakoff effect, 
where 77's are produced by the interaction of a real photon with a virtual photon in the 
Coulomb field of the nucleus. In 1974 at Cornell a measurement based on the Primakoff 
effect gave T{ri — ?- 77) = 0.324 it 0.046 KeV, a value 4(T away from the QED results [43]. 
Recently, a reanalysis of the Primakoff experiment, with a different modelling of the nuclear 
background, brought the value of the width in line with direct measurements, precisely to 
T{r] — ^ 77) = 0.476 ±0.062 KeV [44]. Extraction of the Primakoff amplitude from the data 



-7- 



Figure 1. (a) SOZI and (h) DOZI diagrams contributing to ip ^ PV decays 



is very delicate, because this amplitude interferes with hadronic amplitudes due to vector 
meson [p and w) and axial vector meson bi exchanges. Increasing the energy may help, 
since at very high energies the growth of the Coulomb peak must dominate over the Regge 
behavior of the strong amplitude. After more than 30 years from the Cornell experiment, 
a new experiment to measure the r(?7 — t- 77) decay width via the Primakoff effect has been 
proposed and approved at Jefferson Laboratory, using a 11.5 GeV tagged photon beam 
on two light targets, proton and ^He [45]. The targets have been chosen with the aim 
of minimizing the nuclear incoherent background and enabling a good separation of the 
Primakoff production mechanism from the nuclear coherent background. They estimate to 
reach a 3% accuracy in the measurement of the r/ width, that would yield less than 1° of 
uncertainty on the r/ — ?/ mixing angle. 

As it is well known, all rj meson possible strong decays are forbidden in lowest order 
by C, CP invariance and G-parity conservation. First order electromagnetic rj decays are 
forbidden as well, or occur at a suppressed rate because of involving an anomaly. The first 
allowed decay is therefore the second-order electromagnetic transition rj — )• 77. The decay 
r/ — 7- 37r violates isospin symmetry and it is mainly due to the isospin breaking part of the 
QCD Lagrangian, since contributions from the electromagnetic interaction are strongly 
suppressed by chiral symmetry [46]. The main interest of this decay resides in the fact 
that, in principle, it offers a way to determine the mass difference of the up-down quarks. 
The absolute value of the partial decay width for 77 — )• 37r is experimentally obtained via 
normalization to — >• 77; therefore, a change in one decay width has influence on the other 
[30]. 

A few comments are in order for the analysis of -0 — )■ PV . It was pioneered by Mark 
III in 1985, when they inferred from their data Z"^, = 0.35 it 0.18 [47]. They assumed 
that such decays proceed via singly disconnected diagrams (SOZI) with their strong quark 
line correlations and ignored doubly disconnected diagrams (DOZI). In Fig. l.(a) and in 
Fig. l.(b) we show examples of SOZI and DOZI diagrams. Motivated by the measurement 
oi ij: ^ ^(^(j), which showed the relevance of DOZI-suppressed processes in -0 decays, they 
performed a new analysis [48], including DOZI contributions and any additional component 
as gluonium or radial excitation. The new analysis did not show evidence for non-qq 
components in the r] and rj' wave functions. 

In 2007 Thomas [23] - following the approach of Seiden et al. [49] - investigated the 
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strong ijj —7- PV transitions; he concluded that DOZI contributions are significant, and that 
any gluonium components should play a role similar to that of DOZI contributions. From 
such an analysis he finds that the fit favors a small gluonic component in the r/', with no 
great significance. Without form factors, Thomas finds (f)p = (45±4)° and (f)G = (33 it 13)° 
(i.e. = (0.30 ± 0.21)), whereas with form factors 0p = (46^^)° and 0g = (44 ± 9)°, 
(i.e. Z'^i = (0.48 lb 0.16)). Another phenomenological analysis of ^ — )• PV , without form 
factors, by Escribano [27], finds <j)p = (40.7 it 2.3)° in the hypothesis of no gluonium and, 
allowing for it, (t)p = (44.6 ± 4.4)° with Z'^, = (0.29l[];^^). 

The remaining decays of the list (3.1), (3.3) and (3.4) are charmonium decays into 
^ri^'\ BESII data have better precision than previous measurements; according to the 
hypothesis of no gluonic contribution, SU(3) flavor symmetry and exact OZI rule, they 
extract an angle in the octet-singlet scheme. Their value, translated in the flavor scheme 
according to the relation (2.3), reads 4)p = (32.62 it 0.81)° [50], a quite low value compared 
to other determinations. The extraction of the mixing angle in [50] has been performed in 
a very symmetric -and therefore simplifled-scheme; we observe that just by introducing a 
dependence on a strange/nonstrange factor, the author in [23] finds for the same processes 
and PDG averaged data (including BESII results) values of the mixing angle in line with 
determinations from other processes. If there is any charmonium component in the rj^'\ 
we expect the decays of ip and ^' into 77/^') to be dominated by the magnetic dipole 
transition of charmonium. In that case, it is possible to estimate that the amplitudes of 
the charmonium components of the 77^'^ are negligible, being less that 5% [23]. 
More recent measurements of ^t]^'^ branching fractions have been reported by CLEO-c [51]. 
The last update of the ip — )• 777' branching fraction has been given by BESIII [52] and reads 
B{Tp 777^) = (4.84 ± 0.03{stat) ± 0.24(sys)) x 10-^ which is consistent with the BESII 
value within 1.5a and with the CLEO value within lAa. The ip' — >• 77/'-'^ decays have also 
been observed by BESIII [53], but no new mixing angle estimate has been reported by the 
Collaboration. As far as T(IS') — )■ 777^'^ is concerned, only upper limits are available for 
the branching ratios from CLEO III [54]. 

Since all extractions of the mixing angle involve some nontrivial theory assumptions, 
it is not totally surprising to find different compositions of the wave functions, yet it is 
still frustrating. The best short- or midterm prospects for improvement lie in obtaining 
constraints from more data of even greater variety. 

Let us now provide some estimates of how much future data can reduce most of the 
uncertainties discussed here. 

3.2 Improving the Constraints of the rj — rj' Wave Functions 

The determination of mixing angles and gluonium content is based on measurements. The 
significance of such constraints depends on the experimental uncertainties. Therefore we 
analyze which experimental inputs will best improve our knowledge of the rj — rj' wave 
functions. We start with the PDG 2010 values [55]: 

• The stated (j) ~^ v'l partial width is mainly due to the KLOE measurement in [24]; 
the error is dominated by systematics due to the secondary 77' branching ratio. The 
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(/) — 7- 7/7 branching ratio has been accurately measured by CMD-2 and SND [55]. 

• The r( — )• uj^ partial width of (0.0053 it 0.0005) MeV with a relative error of 9% comes 
from the overall PDG 2010 fit. The relevant experiment has been performed in 1977 
and was based on 68 events [56]. The KLOE-2 Collaboration [41] could measure the 
branching ratio B{rl — t- a;7) more accurately by collecting at least 20 fb~^ of data; 
the limiting factor then comes from the uncertainty in the total r( width, F^/, since 
it is the partial width that matters. 

• The r( fry partial width inferred from the PDG 2010 fit is (0.0568 ± 0.0030) MeV; 
the absolute branching ratio measurement was performed in 1969 by Rittenberg [57] 
based on 298 events. The PDG fit value is slightly lower than the directly measured 
one. Again the error is dominated by the uncertainty in T^/. 

• The latest values on /? — 777 and uj ^ rj^ partial widths are obtained in [58], based 
on SND data on e"'"e~ — )• rj'j: their accuracy is quite comparable to that of the PDG 
2010 fit values. 

In Table 3 we sketch different experimental scenarios. Starting from the present status 
as given by PDG 2010 we analyze the impact various conceivable improvements in the 
experimental constraints would have on the determination of the mixing angle (pp and the 
size of Z^i, the gluonic component in the r]' wave function. We have chosen the radiative 
processes that are common to analyses [23-25] discussed in Sect. 3.1. 

In column I we list the uncertainties in the experimental input values as stated in 
PDG 2010. In column II we indicate the improvement that could be achieved by studying 
rj' — )• ujj with a sample of 20 fb~^ of e"'"e~ — >• cf) events, that KLOE-2 anticipates to acquire 
in the next few years [41]. We assume a selection efficiency of order 20% in the analysis 
of </) — 7- 7/7 with 7/ — >■ and neglect background subtraction. We observe that the 
limiting factor is provided from the uncertainty in the total rj' width. In column III we 
indicate the improvement that could be achieved by reducing the uncertainty on rj' — )■ p7 
of one-half respect to the present scenario; such improvement is also possible after a few 
years of running of KLOE-2 [41]. In column IV and V we indicate the sensitivity to an 
improvement in the determination of the partial widths for (p — )• rj^'^j and for all the partial 
widths, respectively. Among possible secondary decays of </> — )■ 7/7, there are both decays 
T]' — )• pj and ij' — )■ LO'j, whose errors are dominated by the uncertainty on F^/. However, 
the former is more convenient to measure, e.g. at KLOE, since it has a branching ratio 
of almost an order of magnitude larger; also the total p decay width Fp is much larger, 
partially including and obscuring, from an experimental point of view, the total cj decay 
width Fi^. 

Since the partial widths of processes containing rj' and the total width F^/ are corre- 
lated, in column VI we evaluate the impact of the reduction of the uncertainty on F,;^'. We 
assume a future F^/ measurement with 1.4% uncertainty, which is within the possibility of 
KLOE-2 [26]. Such a measurement allows a determination of a nonzero gluonium content 
at 5a, as shown in column VI. The crucial quantities to consider are not the central values 
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for (j)p and Z^,, since they are likely to shift, but their uncertainties. We conclude it is 
most important to reduce the uncertainty in the partial width for r/' — ?• p7; i.e., one has to 
measure both B{r]' — )■ pj) and T^/ more accurately. 



Processes 




(^r/r),, 


{5T/T)ni 




{5T/T)v 




r?'7 


3.5% 


3.5% 


3.5% 


1.7% 


1.7% 


1% 




2% 


2% 


2% 


1% 


1% 


2% 


rj — !• a;7 


9% 


4.5% 


9% 


9% 


4.5% 


1.7% 


f] p-y 


5% 


5% 


2.5% 


5% 


2.5% 


1.7% 




7% 


7% 


7% 


7% 


3.4% 


7% 


cj — >■ 777 


9% 


9% 


9% 


9% 


4.5% 


9% 


4>p 


(40.6 ±0.9)° 


{m.i+Jiir 


(40.7 ±0.7)° 


(40.6t°:=e)° 


(40.4 ± 0.5)° 


(40.1 ±0.3)° 




(0.09 ±0.05) 


(0.13 ±0.05) 


(0.08 ± 0.04) 


(0.09 ±0.03) 


(0.10 ±0.03) 


(0.13 ±0.02) 



Table 3. /; widths from PDG 2010 fits; II: errors on rj' — > lo^ reduced; III: errors on r/ — > fr/ 
reduced; IV: errors on (j) ^ 77^'^ 7 reduced; V: reducing the uncertainties for all the partial widths; 
VI all recalculated in the hypothesis of 1.4% for the rj' full width. 



The situation concerning the rj' full width is somewhat curious at present: PDG 2010 
lists as its best value T^i = (0.194 it 0.009) MeV - with the error including a scale factor of 
1.2 - resulting from an overall fit. Direct measurements from 1979 [59] and 1996 [60] on the 
other hand yield the average T^/ = (0.30±0.09) MeV, which would lead to <i)p = (42. 7;^;°)° 
and Z"^, = (O.OOibO.13). Recently a new measurement has been performed at the COSY-11 
facility: T^/ = 0.226 ±0.017(stat)±0.014(s?/st) MeV; the value of the width was established 
directly from the measurement of the mass distribution of the rj' meson, determined with 
a very high resolution [61]. The present average world value (2011 PDG partial update) 
contains this last measurement and gives F^/ = (0.199 it 0.009); in the global fit to the rj' 
partial widths the correlations among the partial widths do not change significantly. 

Let us observe that the total width P^/ extracted by PDG and the value of the partial 
width T{rj' — t- 77) are strongly correlated, which may create difficulties when the total and 
the partial width are used at the same time, as in the present case of the mixing angle 
extraction. Moreover, the branching ratios of the rj' meson decay channels are generally 
known with a relative precision of more than an order of magnitude better than the present 
accuracy with which P^/ is extracted. 

4 Weak Decays of Charm and Beauty Hadrons 

After many years of strenuous efforts to obtain the rj and rj' wave functions with nontrivial 
bounds why should one not declare "victory" and go on to something else? There are three 
reasons: 

• Professional pride - not to be belittled in Italy and Bavaria. 

• Lattice QCD simulations have just entered the adult period. 
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Figure 2. Spectator diagrams for D^^-j — > rj^'Hv and _B+ — >■ rj'^'Hv decays. 

• Yet there is the most topical reason, namely that knowing reliably the r/ and 77' wave 
functions are an important input for our understanding of several weak decays of 
beauty and charm hadrons. Most crucially we need it for predicting CP asymmetries 
involving i] and r]' in the final states and to understand whether a deviation from 
standard model (SM) predictions can be seen as a signal of physics beyond the SM 
[62-65]. 

The SuperB and Super KEK B factories approved in Italy and in Japan, respectively, 
will produce crucial statistics needed for Sj-g) — rj/vj' X and — )■ rj/rj'X. There is 
a good chance that LHCb will likewise and much sooner. 

4.1 Semileptonic Modes 

Since one expects semileptonic transitions to be driven by SM dynamics only (or at least to 
a high degree of accuracy), their detailed studies teach us lessons on how nonperturbative 
hadronization transforms quark level transitions. We will analyze here what semileptonic 
D and B decays can tell us about the 77 and rj' wave functions and maybe more importantly, 
how our knowledge of those can help us to better understand the decay mechanisms. 

Before going into a more detailed discussion, a few general points should be mentioned. 
The transitions — )• ij^'h^i', — )• •q^'^l'^v and — t- •q^'^l'^v proceed on greatly different 
time scales, since they are driven by weak interactions on the Cabibbo-allowed, Cabibbo- 
suppressed and Kobayashi-Maskawa-suppressed levels, respectively. Yet they can provide 
us with highly complementary information in the sense that they produce the rf^'^ via their 
ss, dd and uu components, respectively. In addition, as explained below, r/^'^ could be 
excited via a, gg component. 

4.1.1 rj^'hu 

According to the heavy quark expansion the so-called spectator diagrams (see e.g. Fig. 2) 
provide the leading contribution to semileptonic as well as nonleptonic charm decays [66]. 

Data on semileptonic decays need to improve greatly before they can constrain the 
physics related to the mixing with the gluonic component. In 1995 CLEO extracted the 
branching fraction B{D^ — )• rf'^'^e^u) from ratios to hadronic decays of the [67]. In 2009 
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CLEO-c presented the first absolute measurement of the branching fraction of B{D'^ 
r/We+z/) [68]; the ratio 



B{Dt r]e+v) 



= 0.36 ±0.14 (4.1) 

CLEO-c 



is in agreement with the previous CLEO results [67]. In semileptonic Dg decays the final 
state hadron has to be produced off an ss configuration; if r] and r/' are pure qq states, i.e. 
Z"^, = = 0, then one finds in the quark flavor basis 



Rd cot'' (j) (4.2) 



with the quantity Rd given by the relative phase space and the ratio of the i] and r/' form 
factors integrated over the appropriate range in g^. To calculate the explicit form oi Rd 
one has to model the q^ dependence of the form factors, but the factorization of the mixing 
angle dependence can help to devise tests of the mixing angle itself (see e.g. [69]). From 
the previous CLEO results [67], using rj and rj' as pure qq states and a pole ansatz for the 
form factors Feldmann, Kroll and Stech inferred 4>p = (41.3 ± 5.3)° [14]; it agrees even 
better than one might have expected with the values given above as extracted from weak 
and electromagnetic transitions. Their value is consistent with the new CLEO data within 
the errors. 

Gronau and Rosner in a very recent paper [70] gave a similar number for T{Df — )• 
r]'l~^L')/T(Df — >• r]l~^v) (among other predictions) applying a very simple model, where Rd 
is inferred from kinematic factors in the quark level; again, r/ and r]' are described as pure 
qq states. 

The transition form factors encode complex hadronic dynamics and momentum de- 
pendence: in [71] they have been expressed through the light-cone wave functions of the 
initial and final mesons. An allowed range for Z^/Z^, is given; at the point = 0, the 
angle (t)p is estimated to be cj)p = (37.7 it 2.6)° and the simple factorized relation holds [71] 



T{Dt r/e+z/) 



Rd cot 4>pcos^4>G (4.3) 



where (pa has been defined in Eq. (2.7). In [71] the value Rd = 0.28 is estimated by 
neglecting the nontrivial dependence on the constituent quark transition form factor, that 
is a conventional approximation in literature, while Rd = 0.23 is estimated by assuming 
a simple monopole dependence. We observe that the mixing angle extracted from (4.3) 
is strongly dependent on the value of Rd] in order to provide a rough estimation of the 
theoretical error we consider an averaged Rd, that is Rd = 0.255 it 0.050. By using the 
experimental ratio of branching fractions (4.1), we estimate Z^, = 0.16 it 0.33exp i 0.23th, 
that is (j)G = (23.3ib25.8ea;pil8.0t/i)°, where the theoretical error refers to the errors on Rd 
and (j)p added in quadrature. The experimental error dominates over the rough estimate 
of the theoretical error and it prevents any conclusion on the gluonic content of the rj — rj' 
system. 
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For the Cabibbo suppressed transitions one finds in the same framework: 

r.)n+^ + ( = tan2 0p 4.4 

In 2008 CLEO-c reported its first measurement of T{D^ — )■ -qe^v) and an upper 
bound on V{D^ — )■ r]'e'^v) [72]. Two years later, the same collaboration presented the first 
observation of D'^ — )• rj'e^v, with branching fraction B{D^ — )• rj'e^u) = (2.16±0.53±0.07) 
X 10"^, and an improved B{D+ rje^v) = (11.4 ± 0.9 ± 0.4) x 10"^ [73]. By using the 
above data and the reasonable assumption Rp) ~ Ro, we estimate from Eq. (4.4) the value 

(^P = (41±4e,p±3t/^)°. 

By including a nonzero gluon contribution, we can parametrize the D"*" ratio as in 
(4.3). However, with the available recent data, the estimate of the angle (j)p can be made 
independently of (pc by taking the ratio 

r(D+ ^ r/'e+i.)/r(Z5+ ^ r?e+i/) " "^"^ ^ 

The left side is given by the recent experimental data quoted before, and we get (pp = 
(40 ±3)°. 

Yet this is not the final word on the experimental or theoretical side. A few years down 
the line we can expect BESIII to obtain an even larger sample allowing a more accurate 
measurement with errors on the angle (j)p going down to about 2%. 

The theoretical situation is more complex. While the spectator diagram generates 
the leading contribution, for a precision study we cannot ignore non-leading ones. The 
so-called "weak annihilation" (WA) process contributes even to semileptonic meson decays 
[66, 74], as can be illustrated most directly for and Dg] see Fig. 3. An analysis 
based on inclusive semileptonic D decays, which considers both the widths and the lepton 
energy moments, shows no clear evidence of WA effects [75]. While WA might affect the 
corresponding inclusive semileptonic width only moderately, it should impact the exclusive 
channels Df — )■ 'q'l'^v and — )• rj'l^u on the Cabibbo-favoured and suppressed levels via 
the r/"s gluonic component. The strength of the effect depends on two factors, namely, the 
size of the gg component in the rj' wave function and on how much gg radiation one can 
expect in semileptonic Df, and decays. Lastly, since the main effect might come 
from the interference with the spectator amplitude, it can a priori enhance or reduce those 
rates. Simple relations such as (4.2) do not necessarily hold any longer. 

4.1.2 B+^r]^'hu 

In — >■ rj^'Hu decays one encounters a situation analogous to that for — >• rj^'Hu except 
that their rates are suppressed by [Kift/Kftl^ rather than |Fcd/Vcs|^ and that the range in q'^ 
is much larger. In passing we just want to mention that one needs to understand their rates 
to determine from T{B — > Xulv) /^{B — t- Xcly) with the hoped-for accuracy of 

about 5% [76]. 

In the spectator ansatz one finds using the quark-fiavor basis 
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D+, D+, B 




n. n 



Figure 3. Valence quarks c/s/d (as well as h/u) emitting two gluons which generate r)/r]' via the 
gluonic component of the wave functions. 



with the factor Rb again describing the relative phase space (much more abundant than 
for D mesons) and the ratio of the integrated form factors. The semi-leptonic form factors 
B — )• rj^'^ have been calculated in [77] from QCD sum rules on the light-cone, to next-to- 
leading in QCD. In frameworks based on QCD factorization the mesons Fock-state wave 
functions enter in the form of light-cone distribution amplitudes. Equation (4.6) keeps 
robust under the dynamical assumptions in [77]. Data on the ratio (4.6) have started to 
appear since a few years. The errors are still quite large, comparable in percentage to the 
ones analyzed in the previous section, and prevent definite conclusions on the glue mixing 
to be drawn. 

In 2007 CLEO has found first evidence for — t- •q'l'^v decay, with branching fraction 
B{B+ r]'l+u) = (2.66 ± 0.80 ± 0.56) x 10"^. This year, also the BABAR Collaboration 
measured for the first time B{B+ r]'l+v) = (0.24 ± O.OS^tat ± COS^j^^t) x lO""^ [78], 
superseding the 2008 upper limit [79]. The BABAR value has a significance of 3. Oct and it 
is an order of magnitude smaller than the CLEO result. 

The same 2007 CLEO analysis also reported a new value of the branching fraction B{B^ — )• 
rjl+u) = (0.44 ±0.23 ±0.11) x lO^'' [80], improving previous 2003 values [81]. The result is 
similar to the newest one by BABAR: B{B+ r]l+u) = (0.36 ± 0.05stat ± OMsyst) x 10"^ 
[78]. By using BABAR data [78], the ratio (4.6) reads 



B{B+ il+u) 



B{B+ r]l+v) 



0.67 ±0.24,f„i± 0.1 l,y,t 



(4.7) 



BaBar 



It is evident that the experimental situation is not yet satisfying, although the previous 
value does not exclude a large gluonic singlet contribution to the ?/ form factor. 

The corresponding ratio involving the Bg mesons, that is B{Bs — )• ri'l^l~)/B{Bs — ?• 
r]l^l~), is also potentially informative on the rj^'^ gluonic content, although experimentally 
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Figure 4. jy*-'^ produced via its gluonic component with the gluons being radiated off different quark 
lines. The transition b ^ s is penguin mediated. 

much more challenging. The results for the branching fractions of modes with two charged 
leptons in the standard model are of order 10"^ — 10~^ [82], suggesting that they are within 
the reach of SuperB and Super KEK B factories. 

4.2 Non-leptonic D and Charmless B Decays 

Although estimates of the mixing angles may come from b ^ c dominated processes, such 
as Bg — )• J/il^r]^'^ (see e.g. [23, 69, 83]), within the SM many charmless nonleptonic B 
decays receive significant or even leading contributions from loop processes, which rep- 
resent quantum corrections. Thus they provide fertile hunting grounds for new physics, 
in particular in their CP asymmetries. Yet to make a convincing case that an observed 
CP asymmetry is such that it could not be generated by SM forces alone, one has to 
be able to evaluate hadronic matrix elements. Such an undertaking is greatly helped by 
knowing the wave functions of the relevant particles. 

Modes such as B ^ rj^'^K and Bs — )■ ry^'^c^ seem particularly well suited in this respect. 
It should be noted that the branching ratio observed for B — ^ rj'K exceeds the original 
predictions considerably for reasons that have not been established yet. Those predictions 
had been based (among other assumptions) on identifying rj' as a pure qq state. Allowing 
for a gluonic component opens the door to diverse decay mechanisms. For example, Kou 
and Sanda [84] suggested producing the rj' meson via its gluonic component with the gluons 
being radiated off different quark lines, see Fig. 4. Having the gluon radiation being emitted 
from a single quark line might be a more favorable dynamical scenario (see e.g. [85, 86]). 

Recent branching ratios values are B{B^ — ^ K^ri) = (1.1 it 0.4) x 10^^ and B{B^ — > 
K'^n') = (6.6 ± 0.4) X 10"^ [55]. The B rj^'^K decays may proceed through tree diagrams 
b — )• uus, but such contributions are colour and Cabibbo-Kobayashi-Maskawa suppressed, 
and by one-loop 6—7-5 penguins. Although the same basic penguin mechanism is expected 
to drive both B — t- rj^'^K and B — t- ttK, the rate of the former is measured to be much 
larger. A possible distinctive contribution are flavor singlet amplitudes that are not allowed, 
if the final state contains only flavor nonsinglet states such as pions and kaons. In flavor 
singlet penguins two gluons couple to the r]' violating the OZI rule and the amplitude can 
get contributions from the pure gluonic component of the r]' . 
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The cases where two gluons are emitted by a single hne (b — )• s g g) together with 
spectator scattering and singlet weak annihilation have been explored in the context of 
QCD factorization (QCDF) [87]. In this approach the constructive interference between 
nonflavor singlet penguins seems already sufficient to enhance the B — )■ rj'K branching 
ratio, without the recourse to flavor singlet contributions; however, due to large hadronic 
uncertainties, a sizable gluonic contribution (up to 40%) to the B ^ rj' form factor cannot 
be excluded. 

In the perturbative QCD approach the impact of the gluonic component on the branch- 
ing ratio - potentially important since it increases the branching ratios B — >• rj'K, while 
decreasing the B ^ K rj one - has been estimated to be numerically very small [88]. The 
phenomenological importance of the r]' gluonic content was instead emphasized in the con- 
text of soft collinear effective theory (SCET) [89]. 

Let us note that the previous exclusive analyses have been performed not later than 
2006, when relevant new data, such as semileptonic B ^ rjlv branching ratios, were not yet 
available. In semileptonic decays there is no enhancement in the B decay into r]' mesons. 
The enhancement is also not observed in Df — t- K~^r]' relative to Df — t- K~^7r^. Recent 
data from BABAR for decays into K* [90] favor an opposite pattern with respect to K, 
namely T{B — t- K*rj') < T{B — ?• K^rf). It would be interesting to check the impact of all 
recent experimental values on the different approaches. For instance in [89], the effort to fix 
the size for the gluonic contribution to the B ^ r]' form factor, in a more constrained way 
with respect to [87], partly depends on fitting nonperturbative parameters to experimental 
data. We have to admit that a quantitative comparison with data is hampered by the 
theoretical uncertainties in nonleptonic decays. 

The large measured branching ratio for B^ — )• Ksr]' by the BABAR and Belle Col- 
laborations [91, 92] greatly improves the usefulness of the decay mode for measuring CP 
asymmetry and to produce significant deviation from the SM prediction. The projected 
SuperB and Super KEK B factories will probe highly nontrivial ranges for new dynamics. 

While it is true that the size of the time-dependent CP asymmetry established in B^ — )• 
r]'Ks conforms well with the SM expectation, one cannot count on an intervention of new 
physics being numerically large there. Having a smallish deviation being significant implies 
good theoretical control over the SM prediction, which in turn requires good knowledge of 
the 77' as well as rj wave functions. 

Finding CP asymmetries in D — )■ r/^'V, rj^'^r], rjcj) and interpreting them as signals of 
new dynamics has two experimental and theoretical advantages: 

• The branching ratios are not very small. 

• The SM can produce only very tiny CP asymmetries. Even small asymmetries pro- 
duce clear signatures for new physics, as long one can control systematic uncertainties. 

5 Summary 

In this paper we have described the status of ongoing investigations, starting from a review 
of the knowledge on the r/ and r]' wavefunctions in terms of quark and gluon components 
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as has been inferred mainly from radiative (p and tp decays. The different determinations 
of the r] — r]' mixing are generaUy consistent, but the message concerning the gluon content 
in the ry' remains ambivalent. The semileptonic , Dg and decays can give other 
constraints to check rj' gluonium role. Moreover a sizable gluonium content could help to 
understand the unexpected high value of the branching ratio for B — )■ rj'Ks decay. 

In conclusion: after many and difficult efforts to understand the rj — r]' wave functions 
it might be seen as "smart" to call it a "victory" and move to another problem. We want 
to emphasize that it is a "noble" goal to improve our understanding of non-perturbative 
effects in QCD, in particular when more "allies" from lattice QCD come to the battle line. 
Furthermore, and maybe even more important, it will help significantly to identify the 
footprints of new physics in CP asymmetries in B and D decays. 
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